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Abstract 

We develop a new transformation theory in quantum physics, where the transforma- 
tion operators, defined in the infinite dimensional Hilbert space, have right-unitary 
inverses only. Through several theorems, we discuss the properties of state space 
of such operators. As one application of the right-unitary transformation (RUT), 
we show that using the RUT method, we can solve exactly various interactions of 
many-level atoms with quantized radiation fields, where the energy of atoms can 
be two levels, three levels in A, V and = configurations, and up to higher (> 3) 
levels. These interactions have wide applications in atomic physics, quantum optics 
and quantum electronics. In this paper, we focus on two typical systems: one is a 
two-level generalized Jaynes-Cummings model, where the cavity field varies with the 
external source; the other one is the interaction of three-level atom with quantized 
radiation fields, where the atoms have A-configuration energy levels, and the radia- 
tion fields are one-mode or two-mode cavities. 
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I. Introduction 



Not only an important method, but also an intrinsic description of symmetry for the 
physical systems, transformation theory is always an interesting topic in physics that has 
acquired a lot of studies. The transformations preserving the measurements are called 
physical transformations. We know that classical mechanics is founded on the symplectic 
manifold, where the transformations of momentum and position preserving the symplecto- 
morphism are the physical transformations, which form a group with composition, called the 
symplectic group. These transformations are also called canonical transformations, since 
they leave both the canonical equations of motion and the Poisson brackets invariant. 

Along the jump from classical theory to quantum theory, the Poisson brackets are re- 
placed by the bosonic commutation relations. The linear transformations of bosons those 
leave the bosonic commutation relations invariant were first introduced by Bogoliubov 1 . 
These bosonic transformations were found to form the same symplectic group as that in 
classical mechanics, even though these two kinds of transformations have different physi- 
cal meanings essentially 2 . The quantum canonical transformations of fermions were first 
introduced by Valatin 3 in his study of superconductivity. Since fermions are the purely 
quantum objects without classical correspondence, the Lie group formed by the fermionic 
transformations is no longer the symplectic group, but a sub-Lie group isomorphic with the 
rotation group 4 . The Bogoliubov- Valatin transformations as the linear quantum transfor- 
mations have wide applications, for examples in the BCS theory 3 of superconductivity and 
in the calculation of blackhole radiation 5 . Recently, the supersymmetric transformation that 
mixes the boson fields with the fermion fields in a unified form advances the transformation 
theory to a big stage 6 . However, some interesting physical predictions raised by the super- 
symmetric theory are still awaiting experiment to test. In quantum gauge theory, the BRST 
transformations of the gauge fields and the ghost fields are in the well-known supersym- 
metric forms 7 , which play an important role in the renormalization proof of the Standard 
Model. In sum, the quantum transformation theory is still under development for various 
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purposes. We notice that the above transformations are the unitary transformations. Since 
quantum theory is set upon the Hilbert space, and the duality of the Hilbert space is defined 
through the Hermitian conjugate, it has been recognized that only the unitary transforma- 
tions do not change the Hilbert space. On the other hand, the nonunitary transformations 
can not preserve the Hermitian duality, thus break the realities of physical observables and 
probabilities. That is why the nonunitary transformations are always rejected in quantum 
theory 8 . 

In a previous paper 9 , we introduced an alternative method: right-unitary transformation 
(RUT), to deal with the two-level Jaynes-Cummings model 10 , which is a basis of the fully- 
quantum description of radiation-matter interaction, and widely used in quantum optics, 
quantum electronics, etc. It was defined in Ref. [9] that if an operator U satisfies the 
conditions: l/W = 1, WU ^ 1, it belongs to RUT. In a strict sense, U is a special 
nonunitary operator that has a right-unitary inverse only. In the matrix representation, 
U can only be the matrix in an infinite dimensional space. Such an operator has been 
recognized as the irregular operator by mathematician, thus obtains few studies, let alone 
its utility in physics. However, in Ref. [9], we found that various JC models can be solved 
exactly by the RUT method. This method not only shows its own merits such as simplicity 
and general applicability, but also leads to a deep understanding of the JC models in essence. 
This work further implies that the transformations applicable to quantum systems should 
not be restricted in the unitary range. Instead, some nonunitary transformations which were 
regarded as the irregular objects might have their particular utilities in physics. 

We know that the unitary transformations do not change the measurements of a quan- 
tum system. But the situation is quite different when using the nonunitary transformations. 
In this paper, we attempt to develop a theory of the right-unitary transformation in math- 
ematics, and discuss its applications in physics. 

This paper is arranged as follows. In Section II, through six theorems, we reveal some 
common properties of the right-unitary transformation (RUT), and conclude a general pro- 
cedure on how to apply RUT to quantum systems. Section II is the foundation of the whole 
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work, it includes three subsections. In Subsection 2.1, we reveal that the state space of any 
operator U e RUT is composed of two independent subspaces { l^ ) } and { j^ 1 ) }, where 
{ |\&°) } is called the Kernel of U, which usually has a finite number of basic elements, and 
satisfies the equation: U{ |\I/ ) } = 0. On the other hand, in the subspace { l^ 1 ) }, U acts 
as a unitary operator. Similar to the unitary transformations, all the operator U G RUT 
in the same Hilbert space form a semigroup, which is called right-unitary semigroup. For a 
certain system with Hamiltonian H, through a theorem we show that if the Kernel of U is 
checked to be isomorphic with a subset of the eigenkets of H, U can be used as a unitary 
transformation to the remaining subspace of H without changing the spectrum. In Subsec- 
tion 2.2, we briefly discuss the application of RUT to the nonstationary system. In order 
to construct the RUT conveniently in the Fock space, we employ the recently introduced 
inverses ( b, tf) of boson creation and annihilation operators (a*, a) 11 in Subsection 2.3, 
and show that the well-known quantum phase operators 12 constructed by a and a) (or b and 
&t) form an Abelian subgroup of the right-unitary semigroup. 

Section III attributes to one of the applications of RUT in physics, where we use the RUT 
method to treat the systems of many-level atoms interacting with the quantized radiation 
fields, where the RUT are constructed by the quantum phase operators. This section contains 
four subsections. In Subsection 3.1, we approach such a model that the atoms have two 
energy levels, and the radiation field is designed to be a one-mode cavity that varies with 
time. Through this model we hope to achieve the goal of controlling the effect such as 
atomic inversion of the system via the external source. This model can be regarded as 
the nonstationary JC model, which is found to exhibit such a property that there exist 
a particular relation between the atomic inversion and the energy exchange of the atomic 
system with the external source. Therefore, by measuring the energy exchange, we can 
understand the dependence of quantities such as the atomic inversion, and the mean photon 
number on time (on the external source). 

In recent years, the case of three-level atoms interacting with quantized electromagnetic 
field have obtained extensive studies. The model where the atoms have A-configuration 
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energy level has been applied to a number of different contexts such as the coherent popu- 
lation trapping, laser cooling 13 , lasing without inversion 14 , and electromagnetically induced 
transparency (EIT) 15 . In Subsection 3.2, we use the RUT method to solve exactly the 
atom-radiation interaction, where the atoms have A-configuration energy levels, and the 
radiation field is designed to be a one-mode cavity field. The case of a two-mode cavity is 
treated in Subsection 3.3, which is further divided into two parts: In the first part 3.3.1, 
we consider such a situation that the detunings of the system are enough large, then we are 
able to construct a unitary transformation to separate the upper level from the system, and 
the remaining two low levels are still treated by the RUT method. In the second part 3.3.2, 
we use the RUT method to solve the system exactly. In Subsection 3.4, we briefly show 
how to construct RUT to solve the matter-radiation interactions where the atoms can be 
three levels in V or = configuration, and up to four levels. Section IV is the conclusion. 

In the present paper, we focus on the RUT method, the analyses of the quantities related 
to the interesting phenomena such as lasing without inversion, and electromagnetically in- 
duced transparency (EIT) are not given. Since there are lots of formula and quantities in 
the paper, the meanings of the symbols are independent in each section. 



II. Right-unitary transformation 

2.1. Stationary case 

Consider a Hermitian operator A with a discrete spectrum 

A|Vn> = A n |^n), n = 1,2, •••,£;, (1) 

the eigenkets \tp n ) are orthogonal mutually, {ip m \ip n ) = S mn . Let us choose an operator U 
to transform A into another frame A'. If A' is still a Hermitian quantity and has the same 
spectrum as A, we call U an applicable transformation to A. There raise two questions 
subsequently: 
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How to transform A into A' by using U? 

Since we require that A' still be a Hermitian quantity, the transformation from A to A' 
is usually chosen as A' = UAW . Then 
What is the basic requirement to U? 

This question is indeed the crux of the transformation theory, and has not a complete 
answer yet. Certainly, as we have pointed out in the Introduction that the unitary transfor- 
mation is the applicable transformation. Here we would like to see how far we can approach 
beyond the unitary transformation. As it required that A' should have the same spectrum 
as A, it is proper to choose an extreme case, A = 1, as an invariant in the transformation. 
This will result in such a requirement to U that l/W = I, in order to maintain the case of 
unity: A = A' = I. We know that in the infinite dimensional Hilbert space, UW = I does 
not mean that U is a unitary operator. In fact, there exists the the following nonunitary 
transformation, 

UW = /, 

(2) 

WU = W^I. 

The aim of present work is to understand how this kind of nonunitary transformation works 
in quantum theory. We here call the operator U satisfying Eq. (2) right-unitary transforma- 
tion (RUT), in order to distinguish it from the other nonunitary transformations. We know 
that RUT are the operators in the infinite dimensional space, they can only be applied to 
the system having the same infinite dimensional state space as RUT, i.e., k in Eq. (1) goes 
to infinity. 

Before carrying a theoretical study on RUT, we first give a simple example of it: quantum 
phase operators 12 , 



where a and a) are bosonic annihilation and creation operators, respectively. It is easily 
calculated that 
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(3) 



Fx = 
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FiFl = 7, 

(4) 

^^=7-10X01. 

Which indicate that F 1 belongs to the right-unitary operators. This concrete example of 
RUT will be helpful to understand the following mathematical approach. The properties of 
the quantum phase operators will be carefully discussed Subsection 2.3. 

We now study the properties of RUT through the following several theorems. 

Theorem I For any operator U G RUT, that W = WlJ has a complete set of eigenkets 
{ l^) }, where { } is constituted by two independent subsets: { l^ ) } and { j^ 1 ) }. 
For the set { |\l/ ) }, U acts as an annihilation operator, 

U\V?) = 0, |*°> G { |*°> }; (5) 
For the other set { l^/ 1 ) }, U acts as a unitary operator, namely, 

W\*)) = |*J>, |*)> G { l^ 1 ) }; (6) 

Proof: For any U G RUT, it follows from Eq. (2) that 
U(l - W) = 0. (7) 
Equation (7) is further left-multiplied by U\ we have 

W - W 2 = 0. (8) 

Equation (8) indicates that the Hermitian operator W has eigenvalues and 1 only, and the 
corresponding eigenstates, denoted by { l^ ) } and { j^ 1 ) }, form a complete set, where 
the brackets are used to represent the case of degeneracy if exists. For any G { l^ ) }, 
using Eq. (7), we obtain U(l - W)\^) = U\V?) = 0, which proves Eq. (5). On the other 
hand, for any |^)) G { l^ 1 ) }, we have UW\^) = W]^) = |^)), namely, U acts as a 
unitary operator in the subspace { l^/ 1 ) }. # 

As the start point of our work, the above theorem reveals the general structure of the state 
space of RUT. Applying this theorem to the phase operator 7\, we have: { l^ ) } = { |0) }, 
and { l^ 1 ) } = { |jfe), k = 1,2, •••oo }. 



The following two corollaries are directly associated with Theorem I: 
First, from Eq. (5), we know that 

Corollary I. a. For any operator U G RUT, there exists at least one vector \ip) ^ to 
satisfies the equation: U\ip) = 0. # 

Second, for an arbitrary vector \9), from that 

WU\0) = U^U [|tf°)(tf°|0) + l^ 1 )^ 1 ^)] = (tf 1 ^)^ 1 ), (9) 

(the brackets are omitted here) we obtain: 

Corollary I.b: For U G RUT, WU is a projection operator that maps an arbitrary 
vector into the unitary subspace of U. # 

It is well known that all the unitary operators in the same Hilbert space form a unitary 
group. Similarly, we have the following group theorem for the RUT: 

Theorem II All the operators satisfying Eq. (2) in the same Hilbert space form a 
semigroup. 

Proof: In this proof, the main task is to prove that for two arbitrary operators Ui, U 2 G 
RUT in the same Hilbert space, that U = U 1 U 2 should also belong to RUT. It is evident 
that UW = I. We now prove WU ^ 1. Let W 1 = U\U X , W 2 = U\U 2 . Then W = WU = 
U 2 WiU 2 . The Corollary I. a tells us that for U 2 , there exists at least one ket vector ^ 
to satisfy the equation U 2 \ip) = 0. Which means W\ip) = 0, and W ^ 1. Therefore, U £ 
RUT. The associative law is easily proved. Since all the RUT have not strict inverses, we 
conclude that all these RUT in the same Hilbert space form a semigroup, which is denoted 
asRUSG. # 

The above two theorems reveal some common properties of RUT, which are independent 
of the concrete structure of RUT. In order to study the structure of RUT in detail, we now 
introduce a useful quantity: the Kernel (K) of the operator U e RUT: 

K = {\$?), i = 1,2,---}, for any |<t>°> E K, [7|<t>°> = 0. (10) 

Evidently, the subset {I 1 ? )} in Theorem I and the arbitrary linear combinations of those 
elements in {I 1 ? )} all belong to K. Therefore, K has infinite number of elements. Here we 
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introduce another quantity--Basic Kernel, which includes only the base ket vectors of K, 
namely, these ket vectors are orthogonal with each other. In most cases, the Basic Kernel 
has finite number of elements. We sometimes simply call the Basic Kernel as Kernel. For 
convenience, we usually denote K = {|$ )}, or K = |$°) in the following presentation. 

The Kernel distinguishes clearly the RUT from the unitary transformation, since for any 
unitary transformation, K is an empty set. In a certain sense, K can be taken as a measure 
of the non- unitarity of the operator U G RUT: the fewer elements K has, the more nearly 
unitary U is. The latter discussion will show that K directly determines the applicability of 
RUT to a physical system. 

Since RUT has been associated with a structure of semigroup, we now prove the following 
Theorem about the Kernel of several RUT's product: 

Theorem III For those Ui, U 2 , • • • , U n , • • • G RUSG, the Basic Kernel of operator 

[/(l,2,---,n) = U r U 2 ---U n (11) 
is exactly 

K={mn)}, Ul\&(n-1)), UluU^(n-2)), U\ ■ ■ • E/J E/J|<&°(1))} , (12) 

where |$°(m)) is the Basic Kernel of U m . (The bracket is omitted here) 
Proof: The proof includes three steps. First, let 

W = C/t(i, 2, • • • , n) U(l, 2, • • • , n) = tf n ■ ■ ■ U\ U\W X U 2 U Z •••[/„, (13) 

using Theorem I, we obtain that 

W\$°(n)) = W [ul\<Z>°(n - 1))] = ■ ■ ■ = W \u\- ■ -U\ C/ 2 f |$ (l))] = 0. (14) 

Therefore, all the elements |$°(n)), Ul\$°(n - I)),---, U^- ■ -U\ U\\$ Q (l)) in Eq. (12) 
belong to the Kernel of C/(l, 2, • • • , n). 

The second step, it is easily proved that all the elements in Eq. (12) are orthogonal with 
each other. We omit the proof here. 



The third step, we will show that an arbitrary element belonging to the Kernel of U is 
uniquely determined by the elements in Eq. (12). Supposing an arbitrary vector \a) that 
satisfies W\a) = 0, by Eq. (13), we have 

W 1 U 2 U 3 ---U n \a) = 0. (15) 

We conclude from this equation that there are only two possibale \a) as follows : 

(i) U 2 U 3 ■ ■ ■ U n \a) G |$°(1)). From Theorem I, we obtain that U 3 ---U n \a) = 
cit/] |$°(1)) + c 2 |$°(2)), where q are the parameters commutative with all C/j. By this re- 
sult, we obtain further U 4 ■ ■ ■ U n \a) = cif/|C/||^°(l)) +c 2 f/||^°(2)) + c 3 |^°(3)). Following the 
same analysis, we eventually get to that \a) = c-JJ^ ■■■U\ f/ 2 t |$°(l)) + c 2 Ul ■ ■ ■ f/ 3 t |$°(2)) + 

n 

• • • c n |$°(n)), which evidently belongs to K, where |q| 2 = 1- 

i=i 

(ii) The other possibility is that Ui • • • U n \a) = 0, while U[ + i • • • U n \a) = rf;|$°(/)) ^ 0, 
where 2 < I. Using the same discussion as in (i), we obtain that \a) = diU^ ■ ■ ■ £// +1 |$°(/)) + 
h d n |$°(n)), which still belongs to K in Eq. (12). 

The above analysis leads to a conclusion that the expression (12) uniquely determines 
the Kernel of f/(l, 2, • • • , n). # 

With the help of above several theorems, we now attempt to establish a connection 
between RUT and A in Eq. (1). It seems that the above properties of RUT are independent 
of A, even though RUT is initiated by its application to A. The following Theorem IV 
expresses a sufficient condition, under which we can apply RUT to A: 

Theorem IV If the Kernel of U, K = {|$j)}, is linearly isomorphic with the set 
{|^), I — 1, 2, • • • , s}, which is a subset of the total eigenkets {\ip n ), n — 1, 2, • • • , oo} of 
the operator A in Eq. (1), then U can be taken as a unitary transformation to the remaining 
subspace of A without changing the spectrum. 

Proof: We complete the proof by two steps. First, from the proposition we know that 

s 

li^i) = J2 du\$i), l < s. By which we obtain further 
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W,)=0, Z = l,2,...,s, 

(16) 

E^m) 7^ 0, ra = s + 1, • • • , oo. 
Above equations indicate that the total state space of A can be divided into two parts, 
which are simply denoted as {IV'j)} an d {iV'm)}, respectively. The subspace {JV'j)} can 
be determined by both of the Kernel of U and A, but the remaining subspace {iV'm)} is 
independent of K. 

Second, we now deal with the remaining subspace {\ip m )}- Let \ip' m ) = WU\ip m ). It 
follows from Eq. (2) that U(\ip' m ) - \ip m )) = 0. Equation (14) implies 

WJ = |VU + I> m M>, (17) 
i=i 

where c m \ are parameters to be determined. \ip' m ) is further left-multiplied by a bra {ipi>\, 
I' < s: 

(Mrtu\i/> m ) = (M*Pm) + j2 c mi(M^i), (is) 

1=1 

from Eq. (16) and the orthogonality of the ket vectors {l^n)} indicated in Eq. (1), we 
obtain c m i = in Eq. (17). Then, \ip' m ) = U^U\il} m ) = \tp m ). This result evidently shows 
that in the subspace {|^ m )}, U acts as a unitary operator. 

It is easy to prove that A' = UAW has the same spectrum as A in the subspace (|-0 m )}. 
The new eigenkets of A', \^ m ) = U\ip m ), are also complete, since 

oo s 

E l*m>(*m| = U{1 ~ E = UU^ = 1. (19) 

m=s+l 1=1 

Moreover, by that U^\^ m ) = \ip m ), w e obtain the eigenkets of A. # 

The proposition in above Theorem IV that the Kernel of U should be isomorphic with a 
subspace of A is very strong. How to weaken this proposition is still under investigation. 

Notice that in the proof of Theorem IV, we have employed the eigenstates of A. Therefore, 
one may wonder how to use the RUT method to obtain the spectrum of A, provided we 
do not know these eigenstates at first. To answer this question, we should keep in mind 
that even when using the unitary transformation to solve a problem, there still is not a 
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widely accepted rule on how to construct a unitary transformation, except some well-known 
problems in Ref. [2]. However, for the RUT method, based on Theorem IV, we conclude a 
general way on how to apply RUT to a concrete problem, there are several steps described 
as follows: 

Consider a Hermitian operator such as A in Eq. (1). To solve the eigenvalue equation, we 
first attempt to construct a transformation U in the same Hilbert space as A to simplify it. 
For the purpose, U is constructed to be a right-unitary operator, its Kernel is easily obtained. 
Then we can directly check whether all the elements in the Kernel (or the combinations of 
the those elements in Kernel) are the eigenstates of A. If they are, the Kernel is called 
proper Kernel, which guarantees that we can directly follow Theorem IV to take U as a 
unitary transformation to transform A into A' = UAW, without changing the spectrum. 

We should notice that in the above steps of using RUT, we need not know the eigenkets 
of A at first. The most important step is to construct a RUT with a proper Kernel, as 
Theorem IV required. In practice, the number of the basic elements in Kernel should be as 
few as possible, because: i) with a smaller Kernel, it is easier for us to check whether it is a 
proper Kernel; ii) From Theorem I, we know that in comparison, the smaller the Kernel, 
the more nearly unitary the RUT is. 

In Section III, we will follow the above steps to treat various matter-radiation interac- 
tions. 



2.2. Non-stationary case 

We now consider such a case that the Hamiltonian of a system represented by A varies 
with time. The evolution of the state obeys the Schrodinger equation (% = 1) 

mm) = i^w% i#>)> = iiM- (20) 

We further assume that there exists another variable R in A(t): A(£) = A(R,t), and the 
state space of R is an infinite dimensional Hilbert space. Consider an operator U G RUT 
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in the same Hilbert space as R, and independent of time. U has a finite Kernel: K = 
(l^) , i — 1, 2, • • • , s}, and a unitary subspace {\^f}) , i — s + 1, s + 2, • • • , oo} as indicated 
in Theorem I. In order to apply U to the system A(i), we present the following Theorem 
which is essentially a revision of Theorem IV: 

Theorem V : For [/ e RUT, if all the elements or the linear combinations of those 
elements in the Kernel of U are exactly the eigenstates of A(t) in Eq. (20), i.e., 

A(*)IO = A n (t)|Q, m = l,2,...,s (21) 

where |£ n ) = X) d nm |\l/^), the matrix (d nm ) has an inverse, and t in Eq. (21) is taken as 

m=l 

a constant parameter, then [/ can be applied to the remaining subspace of A(t) without 
changing the Schrodinger equation. 

Proof: Since the matrix (d nm ) has an inverse, two subspaces {|£ n ), n — 1,2, • • • , s} 
and {|^°), n = 1, 2, • • • , s} are isomorphic with each other. Theorem I indicates that there 
is a complete set of ket vectors {|£ n ), ^ = 1, 2, • • • , s; i = s + 1, s + 2, • • • , oo} for 

W = WU, where the ket vectors are assumed to be orthogonal with each other. Then the 
wave function of A(t) can be expressed as 

s oo 

Mt))= E/m(0IU+ E 9n(t)\K), (22) 
m=l n=s+l 

where the time- dependent parameters f m {t) and g n (t) are to be determined. Taking this 
expression (22) into Eq. (20), by Eq. (21) and the orthogonality of those ket vectors, we 
obtain that the parameters f m (t) satisfy the equation 

^fm(t) = -i\ m (t)f m (t), m = 1,2, (23) 
which is solved as 



fm(t) = /m(0) exp 



t 



i \ \ m (t')dt' 







, m = 1,2,- (24) 



In the remaining subspace {|^), i — s + 1, s + 2, • • • , oo}, since U acts as a unitary trans- 
formation, the Schrodinger equation (20) can be turned into 
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A'(t)\my = ^ t \my, (25) 

00 

where A'(t) = UA(t)W, and ^(t))' = £ g t (t)U\V}). In the transformed frame A'(t), the 

Z=s+1 

parameters g„(t) are obtained to satisfy the equation 

Q OO 

QL9nit) = -i Y, 9i(t)«\A'(t)\*iy, n = s + l,s + 2,...,oo, (26) 

Z=s+1 

where l^; 1 )' = The equations (25) and (26) evidently show that U can be applied to 

the remaining subspace of A(t). # 

The condition Eq. (21) is the key point of this theorem, by which the total state space 
of A(t) can be divided into two parts: one part is determined by the Kernel of U ; the other 
part can be treated by U as a unitary transformation. In the above proof, the subspace 
(l^j 1 ), i = s + 1, s + 2, • • • , oo} has been used to expand the wave function. In fact, for a 
concrete problem, we do not know this subspace at first. Instead, we can directly apply U 
to A(t) as shown in Eq. (25) to obtain this subspace, provided the condition Eq. (21) is 
satisfied already. 

The above U G RUT is assumed to be independent of time, since R is an time- 
independent variable in A(t). One can extend Theorem V into the case of many variables 
without difficulty 



2.3. Inverses of bosonic operators and examples of RUT 

As one of the infinite dimensional Hilbert space, the Fock space has been widely used 
in quantum physics. In order to construct RUT conveniently in the Fock space, we employ 
the recently introduced inverses of boson creation and annihilation operators a* and a, 
[a, a)] = I. Mehta et al. u found that a) has a left inverse b, and a has a right inverse b\ 

bet = atf = 1. (27) 

The properties of b, and their applications to squeezed states, and the Mobius transfor- 
mation have obtained detailed studies in Ref. [12]. In Ref. [9], we showed that b and can 
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be formally expressed in terms of a and at, 

b = (±r)a, 6t = a t ( J_). (28) 

From Eq. (28), we obtain further: btf = 1/aa*. Using this relation, we arrive an interesting 
result: a and a* can be formally expressed in terms of b and too, 

We know in the scheme of second quantization, an arbitrary operator can be expanded by a 
and , namely, a and a* form a complete operator set. Equation (29) implies that any oper- 
ator can also be expanded by b and . Therefore, b and form a complete operator set too, 
which is connected to that of a and a} by the nonlinear transformations (28) or (29). 

Equation (28) significantly simplifies the calculation in the representation of b and . 
Various results in Ref. [12] such as the analytic studies of b and in the representation of 
Bargmann space are easily obtained by using Eqs. (28) and (29). For example, using Eq. 
(28), one can prove the following useful formulas: 

a k a) k (N + l)(N + 2)---(N + k) (N + k)V 1 ' 

oo oo k— 1 

tf k a k = Y J tf k \n)(n\a k = ]T \n + k)(n + k\ = 1 - ^ \n){n\, (31) 

n=0 n=0 n=0 

where N = a^a. 

We now look at the well-known operator — Phase operator (phasor) Fi as shown in Eq. 
(3), which is initiated from the quantization of phase factor in quantum mechanics 12 . In 
this paper we will use only the photon-lowering property of these phase operators, and that 
the phase properties of the electromagnetic field are not calculated. The higher order phase 
operators are defined as: 

F k = (F 1 ) k = -j^=a k , Fl ee (Fl) k = a tfe ^^. (32) 

In the representation of b and , F\ and Fl are found to have the same forms as Eq. (3), 

15 



F 1 = -4=6, F\ = 6 t ^L=. (33) 
By Eq. (30), we have 

F k = -r^=b\ Fl = 6 tfe ^==. (34) 
With these preparations, we now prove the following theorem: 

Theorem VI: All the phase operators defined by Eqs. (3) and (32) form a subgroup 
of the right-unitary semigroup (RUSG): {1, F k , k = 1, 2, • • •} C RUSG. 

Proof: It is clear that F k Fl = I. On the other hand, using Eq. (31), we obtain 

k— 1 

Fl F k = I-J2\n)(n\^L (35) 

n=0 

These mean that F fe e RUT. Moreover, = e RUT. We therefore conclude that 
k = 1,2, ■, oo} form an Abelian subgroup of RUSG. The Kernel of the element F k is 
directly obtained as K(F k ) = {|0), |1), • • • \k - 1)}. # 

Using this explicit example of RUT, one can easily check the theorems presented above. 
We are not going to dig deeper the theory of RUT here. It is important to see how to apply 
RUT to the physical system. 



III. Application of RUT to the system of many-level atoms interacting with 

quantized radiation fields 

Under the rotating-wave approximation, the interaction of two-level atoms with the 
quantized radiation field is described by the Jaynes-Cummings (JC) model 10 , which has been 
extensively applied in quantum optics, quantum electronics, etc. Various modifications and 
generalizations to the original JC model have been made to approach quantum effects such 
as quantum collapses and revivals of atomic coherence 16 , squeezing phenomenon 17 and so 
on. These JC models have a common attractive property that they all can be solved exactly. 
Since the super symmetric structure was found to exist the JC model, the JC model is viewed 
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as a generalization of the supersymmetric harmonic oscillator system, and its solvability may 
be interpreted in terms of supersymmetric breaking 18 . Using a deformed oscillator algebra, 
Bonatsos et al. 19 gave a unified solvable formulation of various JC Hamiltonians. Yu et al. 20 
pointed out further that there embeds an unusual su(2) algebraic structure in these JC 
Hamiltonians. 

The method mentioned above are valid to the two-level stationary JC models only. It is 
difficult to extend these methods to the time-dependent JC models, or to the case where the 
atoms have higher (> 2) energy levels. In a previous paper 9 , we employed a right-unitary 
transformation (RUT) to solve exactly a generalized two-level JC model and pointed out 
that there exists a geometric phase in the model. 

In this section, we will show that the RUT method developed above can be generally 
applied to various forms of atom-radiation interactions, where the atoms can have two, three 
and higher energy levels in different configurations, and the radiation field can be one-mode 
or two- mode cavity. This section includes four subsections. In Subsection 3.1, by construct- 
ing a RUT, we follow Theorem V to solve such a model that the two-level atoms interact 
with a radiation field which is dependent on time. This model is called nonstationary JC 
model. The importance of three-level matter-radiation interactions has been pointed out 
in the Introduction. However, some interesting characteristics such as the supersymmetric 
structure, and the su{2) structure are not embedded in the three-level models. There still is 
not a generally accepted method to treat these models, according to our knowledge. In Sub- 
sections 3.2, 3.3 and 3.4, following Theorem IV, we will show that various interactions of 
many-level atoms with one, or two-mode cavity can be unitedly treated by the RUT method, 
where the energy of atoms can be two-level and three-level with A, V and ^-configurations. 
The procedures of the treatment are as simple as those in the two-level JC model. 
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3.1. Two-level nonstationary Jaynes-Cummings model 



The approach of the stationary two-level JC models by using the RUT method can 
be found in Ref. [9]. In this subsection, we consider such a case that the radiation field 
interacting with the two-level atoms varies with the external source. Through this model, 
we hope to achieve the goal of controlling the effects such as atomic collapse and revival, 
and the statistics of photon number through the external source. Under the rotating-wave 
approximation, the Hamiltonian of the system with density p(iV)-dependent multi-photon 
(I) interaction has the following general form (h — 1): 

H(t) = coal a + ^cu a 3 + <J l p(N)>y*(t)<T- + p(N)a l -y(t)a+, (36) 

where ui and u are the field and atomic transition frequencies, respectively. a± = (<7i ± 
ia 2 )/2, where a%, 02, 03 are three Pauli matrices. j(t) represents the change of radiation field 
with the external source, j(t) 7^ 0. The detuning A = uo — leu should satisfy the condition: 
|A| <C uj ,uj, in order to preserve the reliability of the rotating-wave approximation. The 
evolution of the system with time is described by the Schrodinger equation 

H(t)\*(t))= i ^(t)), |¥(0)> = |*o>. (37) 



To solve this equation, we construct the following operator U 

(1, 0) 

U = , (38) 

where F t = ^= b l is the phase operator given by Eq. (34). It should be mentioned that 
this U is expressed in the two- dimensional representation of the Pauli matrices, where the 
parameters are phase operators. From Theorem VI, we know that U has the properties 



UW = I, 

(39) 



U, \ 



l-l 



0, 1- E \k)(k\ j 



which indicate that U belongs to RUT, and the Kernel of U is 





m = 


f ^ 


> 


K = < 




, fc = 0, 1, - • - Z — 1 > 











(40) 



One can check that the vectors |^°), k — 0, 1, • • • I — 1, in K are exactly the eigenstates of 
H(t), where the eigenvalues 



E° k = uk - uj /2. 



(41) 



These results evidently show that U is covered by Theorem V. We now follow Theorem V 
to solve the equation (37). The wave function \^(t)) can be divided into two parts: 



|¥(f)> = |¥°(f)> + \*\t)), 

where \^°(t)) is related to the Kernel and given by Eq. (24) as 

l*°(*)> = E/^ (a,fe_a,0/2)t l^)- 



(42) 



(43) 



k=0 



the other part | x I /1 (t)) is determined by the unitary subspace of U. To obtain |^/ 1 (t)), we 
make a transformation: H'{t) = UH(t)U\ that is 

(cuN+^o, gi (N) 7 (t) \ 



H\t) 



(44) 



[ gi (N)Y(t), u(N + l)-±u ) 



where gi(N) = p(N)[(N + l)\/N\] 2 . The Hamiltonian now turns out to be a function of the 
photon number N, where the creation and annihilation of photon have been erased by the 
transformation U. 

We know that in the original frame H(t), the transition from the atomic state 2 to the 
state 1 is induced by the dipole term a" p(N) r y* (t) , and the transition is always achieved by 
creating / photons. On the other hand, the transition from the state 1 to state 2 is achieved by 
the annihilating the same / photons. However, in the new frame H'(t), all of these transitions 
are caused not by the creating or annihilating photon. The transition induced by dipole can 
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only happen between those states having the same number of photons, as shown in Fig. 1. 
Therefore, the states with different photon numbers can be treated independently. 
H'(t) is further rearranged into the following compact form, 

H'(t) =B -f2-M(t), (45) 

where B = lj(N + 1/2), B 1 (t) = Re[y(t)]gi(N), B 2 (t) = -hn[y(t)]gi(N), B 3 = (lu - l)/2, 

— » — » 

jl = —a. 

The above result [Eq. (45)] reveals that in the transformed frame, the Hamiltonian H'(t) 
of the nonstationary JC model describes the motion of a spin-| atom under an effective 
magnetic field B(t), which is a function of the photon number N and time t. 

The evolution of subspace ^(t))' is described by the Schrodinger equation (23) 

H'm'it))' = £\*\t)y, (46) 



dV 



oo 



where ^(t))' = £ 

fe=0 



( a k (t)\k)^ 



the coefficients a k {t) and Pk(t) satisfy the equations 



y Pk(t)\k) j 

dk(t) = -ia k (t)(uk + cj /2) - ip k (t)gi(k)-f(t) : 
$k(t) = -tP k {t)[uj(k + l) + u /2] - ia k (t) gi (k)r(t), 

which are further arranged into 

d k {t) + d k (t)A' k [l, 7 (0] + a k (t)B' k [l, 7 (t)] = 0, 
Pk(t) +$ k (t)A' k [l, Y(t)]+P k (t)B' k [l, 7 *(t)] = 0, 



(47) 



(48) 



where A' k [l, 7 (f)] = iu(2k + l)-^(t)Mt), B' k [l, 7 (t)] = fl?(*)l7(*)l 2 -»(wfc + w /2)7(t)/7(*) - 
to 2 (k 2 + Ik) - ujuq/2 + ul/A. 

For a given j(t), by solving the above two equations, one can obtain the evolution of the 
system with time. For instance, we choose a special case, 
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lit) = e M , 



(49) 



to approach the dependence of some physical properties in this model on A. With this 
particular j(t), a k (t) and f3k(t) are solved to be 



where 

= -\[u{2k + I) - A] ± Jlfa-ul + Xy+lgtk)]*, 



(50) 



CO 



± 



(51) 



/i fc = (wfc + «f + ^)/ft(A;), 



and the coefficients and (3 k are determined by the initial state in the transformed frame: 

(52) 



fc=0 



(a£ + a k )\k) 
{(~h k at + K 1 a k )\k) 



If the initial state in the transformed frame is chosen as 



then the coefficients a k and fx k , v k are related with each other by the equations 



1^(0))' = £ 

k=0 



V 



(53) 



/ 



a 



+ _ l 



{(J>k — h k u k ), 



(54) 



For convenience, we still use the initial state as Eq. (52) in the later calculations. 



Since the initial state is arbitrarily chosen, the equation (50) is in fact the general solution 
of the non- stationary JC model in the transformed frame. From Theorem V, we know that 
the solution can be directly turned into the untrans formed frame by 



(55) 
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Combining this solution and the Kernel part in Eq. (43), we obtain the exact solution of 
the non- stationary JC model for an arbitrary initial state as 



?1=0 



( cO 


( 

oo 


+ E 

k=0 

\ 


\\n)j 



.,+ „ibl, t 



iui, t 



\k), 



• (56) 



Using the normalization condition: (\l/(t)|\E'(t)) = 1, we know that all the coefficients should 
satisfy the equation 

l— 1 oo 



ElM 2 + E[(i + ^)l« fe + l 2 + (i + ^- 2 )l« fc - 

n=0 fc=0 



= 1. 



(57) 



The above solution (56) shows that the Kernel term is independent of "f(t), thereby a 
trivial term. Our interest is in the second term in the solution (56), because this term is not 
only related to j(t), but also determined by the unitary subspace of U. We now evaluate the 
following physical quantities associated with the non-stationary JC model by using RUT. 



i) Statistics of the photon number 

Let n^t) = (*(i)|(<7 ± a 3 )N/2\^(t)) , where a is a 2 x 2 unit matrix, a 3 is the third 
Pauli matrix. For n + , since U [(a + a 3 )N/2] = (<r + a 3 )N/2, we can treat it by RUT in 
the transformed frame, 



n + (t)=[(* 1 (t)\rf] (a + a 3 )N/2 [u\*\t))] = £ k[\at\ 2 + |«^| 2 + &(f) 
where 



k=0 



(58) 



&(f) = ata~ k e-*K + a+ar/VK 



(59) 



By the same way, we have 



i— 1 oo 

n~{t) = Y, n \fn\ 2 + E(* + [Kl 2/ ^ + KIV - 

n=0 fc=0 



(60) 
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It is proper to refer to n + (t) and n~ (t) as the mean photon number of the field for the cases 
where the atom is found in the excited and ground states, respectively. is an important 
quantity in this model, it measures the collapses and revivals of the atomic inversion in a 
single mode with the frequency 

A, = u+ - uZ = ^o-cu/ + A) 2 + %(£;)] 2 , (61) 

Afc is usually called Rabi frequency, where the parameter A obviously affects the Rabi fre- 
quency. 



ii) The atomic inversion 

The atomic inversion in the JC model is illustrated by the quantity a 3 = (t)\a 3 \^I (t)) . 
If we take the two-level atoms as the neutro particles with spin-|, then <7 3 /2 measures the 
mean spin value which varies with time. Since 03 is invariant under U, namely, Ua^W = 03, 
we can directly calculate this quantity in the transformed frame. Using the normalization 
condition Eq. (57), we obtain 



k=0 



a 



+ 12 



+ KI 



- 1. 



(62) 



This result shows that the atomic inversion is measured by the same as that appears 
in the statistics of photon number. In order to evaluate explicitly the effect of A in the 
nonstationary JC model, we choose the following special initial state without the Kernel 
part for an example: 



l*(0)> 



00 

k=0 



I 



-M 2 /2 



/2k\ 



\k) 



\ 



V V2k\ 



\k + l) 



(63) 



For this initial state, we have (see Appendix 

(A + A)<7,(*0 



k\ ) (A + A) 2 + 4^ 2 (fc) 



cos 



^(A + A) 2 + Agf(k) t 



(64) 



k - 

x e 



k\ 



1 (A + \) 9l (k) 

2 + (A + A) 2 + 4^ 2 (A;) 



(65) 
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where the detuning A = u — cul, x = \z\ 2 . With these results, cr 3 is obtained to be 



ff 3 = 



00 / x k e~ xS 



k=0 



k\ 



4(A + %Kg 2 



y/(A + \)* + 4g?(k)t 



(66) 



Equation (66) expresses exactly the collapse and revival of the atomic coherence in this 
nonstationary model, where A plays the same role as the detuning A. By changing A, we 
can effectively control the atomic inversion of the system, thus meet our original purpose 
of studying the properties of the atomic system via the external source. If it is designed: 
A —A, the atomic inversion disappears for the initial state (63). On the other hand, even 
if A = 0, from Eq. (66) we know that the atomic inversion still exists for the external source 
A^O, 



^3 = 

k=0 



00 / x k e -x^ 



k\ 



■ sin 



V + 4g?(k) 



-^\* + 4g?(k)t 



(67) 



iii) Energy 

The mean energy of the system is E{t) = (^(t)\H(t)\^(t)). For the state (56), we divide 
E(t) into three terms, 



E(t) = (*°(t)\H(t)\*°(t)) + U^(t)\Hm°(t))+h.c] + (*\t)\H(t)\*\t)), 



(68) 



(69) 



where the first term in E(i) is easily obtained to be 
E°(t) = l -£(con-u /2)\f n \ 2 . 

n=0 

From Theorem I, we know that the Kernel are orthogonal to the unitary subspace, which 
implies that the second term in above E(t) is zero. The third term can be treated by RUT 
as 



E\t) = [(^(t)^] [UH(t)U ] ] [U\<l>\t)) 



(70) 



A simple calculation gives 
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E 1 (t) = Y,{K\ 2 [hl\-(l + hl)u 



k=0 



+ \ai 



\hf - (i + hf)^} -xut)}- 



(71) 



The variation of the mean energy E(t) with time is still dominated by which measure 
the energy exchange of the radiation-matter system with the external source. 

From the above equations (62) and (71), we abstract the following interesting relation 
between the mean energy and the atomic inversion, 



dE 
~dt 



A da 3 
2~dt' 



(72) 



This relation expresses such a distinct property of this model that by measuring the change 
of energy with time, we can obtain the information such as the dependence of the atomic 
inversion and the mean photon number on time. However, this property does not mean that 
the atomic inversion is caused by the external source. In fact, when A = 0, E turns out to 
be independent of time, but the atomic inversion shown in Eq. (67) still exists generally. 



3.2. Three-level atoms interacting with one-mode cavity 

We now concentrate on the three-level atom-radiation systems. The three-level atoms 
are classified by the configurations of their energy levels. Generally, there are three kinds of 
configurations A, V and =, as shown in Fig. 2, where the atoms with A-configuration energy 
level have been widely used in the subjects such as the coherent population trapping, laser 
cooling 13 , lasing without inversion 14 , and electromagnetically induced transparency (EIT) 15 . 

We know that the exact solvability of various two-level JC models is attributed to a 
unified formulation of JC models by a deformed oscillator algebra 19 , or by a su{2) structure 20 . 
However, there is not a unified solvable formulation to these three-level systems. Start 
from this subsection, we will show that by the RUT method, we can solve various three- 
level models mentioned above. To avoid the paper being too long, almost all the physical 
properties associated with the solutions are not discussed in this paper. 
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In this subsection, we consider a system of three-level atoms with energies £1, e 2 and £3 in 
the A-configuration, which interact with a one-mode cavity field as shown in Fig. 2. a. The 
interaction is considered to be density [pi(N), p2(V)]-dependent, and in the multi-photon 
l 2 ) form. Under the rotating- wave approximation, the Hamiltonian of the system is 
written as 



H = Y^EiSu + ucta + W h Pl {N)S 12 + p 2 (N)a h S 23 + h.c], 



(73) 



i=l 



where uj is the frequency of the radiation field, SV,- are the atomic operators given by Sij = 
i, j = 1, 2, 3. We notice that various three-level Hamiltonians with one-mode cavity in 
the literature are covered by this general one. To solve this Hamiltonian, we now introduce 
the following operator U : 



( F h 



U 



V 



(74) 



1 

F l2 

where and Fi 2 are two phase operators. From Theorem VI, we know that the operator 



U belongs to RUT. The Kernel of U is 



K 



i«i> 


V ) 



, mk 2 )) 



( ^ 





V l*2> ) 



ki < h, k 2 < l 2 



> . 



(75) 



Theorem IV tells us that if this Kernel is isomorphic with to a subset of the eigenkets of 
H, then the operator U can be applied to H . It is easily checked that above IV'i(^i)) an d 
IV'"^)) happen to be the eigenkets of H . The eigenvalues corresponding to these eigenkets 



are 



(76) 



E$(k 2 ) =e 3 + u;k 2 , 

where we have assumed that the energy spectrum is nondegenerate between two sets 
{\ipi{ki))} and {\ip 2 {k 2 ))} ■ These results indicate that the operator U is covered by Theorem 
IV, and can be applied to the remaining eigenket set of H besides those in Eq. (75). 
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Let H' = UHW, a direct calculation gives 



where Hi is 



(77) 



#1 



&i ^i(iV) 
9i(N) -(h + h 2 ) g 2 (N) 
£ 2 (W) /i 2 



\ 



(78) 



where hi and h 2 are two constants: /ii = (2e 1 — e 2 — s 3 )/3 + (2/ x — l 2 )u/3, h 2 = (2e 3 — ei — 
e 2 )/3 + (2/ 2 - Zi)w/3; and ^(JV) = pi(iV)[(iV + h)\/N^, g 2 (N) = p 2 (N)[(N + l 2 )\/N\]l 
Notice that Hi is constructed to be a traceless operator matrix, this form is easier to compute 
its eigenvalues. 

The above results show that in the transformed frame, the Hamiltonian of the three-level 
matter-radiation interaction system becomes a 3 x 3 matrix which depends on the photon 
number iV only. For the ket vector with a single photon number, this Hamiltonian simply 
describes the usual three- level stationary system, which can be solved exactly. We note that 
in the transformed frame, the transitions induced by dipole occur only between the states 
containing the same photon number. This picture is similar that in the two-level system. 

Assuming the eigenket of Hi as 



|*(n))' 



\n). 



(79) 



where $ = Pi(n), i — 1,2, 3. The eigenvalue A of Hi is determined by the equation det — 
A | = 0, that is 



A 3 + \p + q = 0, 



where 



(80) 
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p=- [gl(n) + gl(n) + h\ + h 2 2 + h x h 2 \ , 

q = h x gl{n) + h 2 gf(n) + (h x + h 2 )h 1 h 2 . 
The solutions of Eq. (80) are given by 21 , 

Ai(n) = A+ + A-, 

[\ Vi {n) = ±(A + + A-)±&*(A + -AJ), 
where A± = A±{n), 



(81) 



(82) 



Aj 



(83) 



Combining these Aj with the diagonal term in H', we obtain the eigenvalues of the Hamil- 
tonian as 



(84) 



£i(n) = |(ei +e 2 + e 3 ) +u(n+ ^) + Ai, 

E 2 , 3 (n) = |(£i + e 2 + e 3 ) + (n + ^) + A 2 , 3 . 

We assume that the energy spectrum is nondegenerate, then, the eigenkets of H' correspond- 
ing to the above eigenvalues are orthogonal mutually, which are obtained to be 



l*i(«)>' 



9i{n) \ 
Xi-hi 

1 

92(n) 
\ \i-h 2 ) 



\n), i = 1,2,3, 



(85) 



Xi—hi 



ga(n) V 



1/2 



. The above 



where 0( n ) are the normalization factors: C«( n ) - 
solutions of the three-level system with A-configuration is in the transformed frame. From 
Theorem IV, we can directly obtain the eigenkets in the original frame H by the equation: 
\%(n)) = W\%(n))', that is 
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Xi-h-i 



\n + h) 



Ci(n) 



\n) 



1,2,3. 



(86) 



V 



92(n) ' 
\i-hi_ 



In + h 



) 



For an arbitrary initial state \cp(0)), the evolution of the state with time becomes 

h-l 2 oo 3 



k=0 i=l 



-^ {k)t mk)) + EE^(«)^ (ri)i I^W>. 

«=0 j=l 



(87) 



where Q(k) = (^(k)\ip(0)), and Dj(n) = (Vj(n)\<p(0)). 

Now we conclude that the total eigenstate set of the three-level atom-radiation system 
is constituted by two subsets as in Eqs. (75) and (86), where the subset in Eq. (75) is the 
kernel of U, but for the other subset Eq. (86), U acts as a unitary operator. 

We recall that for the two-level JC model, there generally exists a relation H' ~ a ■ B. 
However, for the above three-level system, from the expression of H', we know that the 
system can not be taken as a spin-1 particle interacting with an external magnetic field, 
that is, the relation such as H' ~ S • B does not exist here. Therefore, there is no concept of 
spin in the A-configuration system, and many methods developed for the two-level JC model 
do not fit for the three-level case. However, for the RUT method, one may notice that the 
procedure in solving the three-level model is exactly the same as the procedure presented in 
Ref. [9] in solving the two-level JC model. 

To investigate the dynamics of the above system, we usually use the method of density 
matrix. Supposing at t — the density operator as 



fei,fe 2 =0 



(88) 



Using Eq. (87), one can directly obtain the evolution of the density operator with time 



by the equation p(t) = e p(0)e . With p(t), one can calculate some quantities such 
as the atomic inversion, population trapping in this model. We now choose such a special 
initial state that the atoms are in the level % — 1, and the photon state is represented by the 
coherent state. Then p(0) is written into 
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p(o)= E 



00 e~\ z \ 2 z kl z* k2 



ki)(k 2 \S u . 



fci,fe 2 =o \J(ki\k 2 \) 
Make use of Eq. (87), we obtain 



(89) 



(90) 



where 



_i£J_ t 



00 3 



= E ^Tff e-^<*>*|^(Jb)> + E E 



e 2 z n+h gi(n)e' 



-iEj(n)t 



k=o Vk ] - 



n =oj=i y] (n + h)\Q(n)[\j(n) - h] 



(91) 



Using p(£), we can obtain the probability of finding the atom in the state i = 2 which is 

00 

initially in the state i — 1: Pi^. 2 (0 = Z) — 2, A;|p(t)|/c, i = 2), 



fc=0 



Pi^it) = E 



00 e~ x x k+l1 



* gi{k) e - iA '^ 



(92) 



where x = \z\ 2 . We know that the above transition (1 <-> 2) is induced by dipole. A more 
interesting quantity is Pi_ >3 (t), which is induced by two dipoles: (1 2) and (2 <-» 3), 



^1-3(0 = E 



r (k + h)\ 



E 



gi{k)g 2 {k) e-^ 1 
^f(k)[\ t (k) - h^k) - h 2 ] 



(93) 



The statistics of the photon number at the state i = 3 is obtained to be 



ns(t) = E 



g 32 fc I i 1 



* gi(k)g 2 (k) (k + l 2 ) e-^W* 



(94) 



^<?(*)[Ai(*)-/n][Ai(fc)-/i2] 
The atomic inversion of the system under the special initial state is S 3 (i) = (5*3), which is 
explicitly given by 



* W = ?<* + W . 



(95) 
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3.3. Three-level atoms interacting with a two-mode cavity field 



In this subsection, we consider a system of three-level atoms with energies S\, e 2 and e 3 , 
which interacts with a two-mode cavity field: a pump mode (ai, a\) of frequency u 1 and a 
Stokes mode (02, a 2 ) of frequency uj 2 as shown in Fig. 2. a. The interactions are generalized 
to be multiphoton (Zi, Z2) forms, and dependent on densities pi(AT 1 ), p 2 (N 2 ), respectively 
Therefore, the Hamiltonian is 

H = H + H 1 , (96) 
where 

3 . . 

H = SiSa + u&iax + co 2 a 2 a 2 , 

1=1 (97) 
E x = a\ ll p 1 (N 1 )S 12 + P2 (N 2 )a l iS 23 + h.c. 

This Hamiltonian covers various special cases in the literature. In the above Subsection 3.2 
we have shown that the solution of three-level case is much more complicated than that of 
the two-level case. To avoid the complication, one usually treats H by the perturbation 
theory in the interaction picture. Especially in the case that the detunings of the two modes 
are very large, the upper level can be eliminated adiabatically from the three-level system, 
then the system is reduced into a simple two- level case with two quantized modes 22 , where 
the Stark shift terms appear and give arise to some interesting physical effects 23 . 

We divide this subsection into two parts. In the first part 3.3.1, we consider the case that 
the detunings of the system are enough large, then by constructing a unitary transformation, 
we separate the upper level from the system to the first order approximation, the remaining 
two low levels are still treated by the RUT method. We further discuss the relation between 
the so-called "dressed" states and "bare" states from the viewpoint of transformation. In 
the second part 3.3.2, we use the same RUT method to solve the system exactly. 
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3.3.1. Approximate treatment in the case of large detunings 

In the case that the detunings of two modes are very large, we accordingly introduce the 
following unitary transformation Q, 



Q = exp(X), 



(98) 



where 



X = 



-<W 1 ,1 pi(Wi)S'i2 + 5 2 p 2 (N 2 )a l iS 23 \ - h.c, 



(99) 



where 5\ and 5 2 are two parameters to be determined. Under Q, the Hamiltonian Eq. (96) 
is transformed into H' = QHQ^. Using the Baker-Hausdorff formula, H' becomes 



H' = H + H 1 + [X, H ] + [X, H,} + l -[X, [X, H \\ + ±[X, [X, H,]] + 



(100) 



If we choose Si = A 1 1 = (e 2 — S\ — WiA) 1 , S 2 = A 2 1 = (e 2 — £3 — oo 2 l 2 ) 1 , where A x and 
A 2 are two detunings, then 



[X,H ] = -H,. 



(101) 



For the large A,, % — 1, 2, we expand if' up to the first order of 1/Aj as 
H'^H + [X, Hi] + l -[X, [X, tf ]] = #0 + \[X, H,], 
that is 



(102) 



i=i 



if' = £^ala, + £l - — pKN, - l^a 1 ! 1 
I Zii 



S11 + 



£3-^P2(^2-i 2 )a^a 2 2 
^2 



'33 



+ 



\(^ + -^)Pi(Ni-h)p2(N 2 )a\ h a l i 



e2 + ^(A^af 1 + -^p 2 2 (N 2 )a l i4 h 

l±l IA 2 



Sis + h.c. 



?22- 



(103) 
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We have two comments on the above results: 

i. The Hamiltonian H' indicates that if the atomic states are in the form 



\v>Y 





M 



(104) 



the term about S22 in H' does not contribute. Then the system is reduced into a two- 
level atomic system. The two terms —-^p\{Ni — h)a\ l1 a l i Su and — -^p\{N 2 — h)^ 2 °2 ^33 
appearing in the Hamiltonian are called Stark effect terms, which are in terms of the photon 
numbers. 

ii. We regard the operators in the original frame H as "bare" operators, and the atomic 
states as "bare" states, then the operators in the transformed frame H' are "dressed" oper- 
ators, and the atomic states are "dressed" states. The state Eq. (104) is in fact a "dressed" 
state. We know that all the states prepared in experiment are "bare" states. Therefore, 
turned back to experiment, the "bare" state corresponding to Eq. (104) is 

/ , - \ 



\<p) = e- x \<p)' 



1^3) 



(105) 



This result shows that the atomic state i — 2 still contributes in the original system. 

We now solve the Hamiltonian H' by the RUT method developed above. Consider an 
operator V 



V = 



F,(D 





1 
V ^(2) j 



(106) 



where ^(1) 



-1/2 



a[\ F l2 {2) 



a 2 2 are two phase operators. Thus V 



belongs to RUT, and its Kernel is 
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K = { 



^ \h, k 2 ) ^ 




, \^M,k' 2 )) 








\k[, k' 2i 



ki < li, k' 2 < l 2 



(107) 



We obtain that \ipi(hi, k 2 )) are the eigenkets of H', the eigenvalues 



are 



k 2 ) — Ei + oJih + u 2 k 2 ; 



1^2(^1, k' 2 )) are the eigenkets of H' with the eigenvalues 



(108) 



E°(k[, k 2 ) = e 3 + ojiK + oj 2 k 2 . 



(109) 



Where we have assumed that there is no degeneracy between E^{k\,k 2 ) and E 2 (k' 1: k 2 ). 
These results indicate that V is covered by Theorem IV, and can be applied to H' to obtain 
the remaining subspace besides Eq. (107). We make the following transformation, 



H" = VH'V\ 



that is 



(110) 



H" = B Q (N U N 2 )(S n + S 33 ) + £3(^1, N 2 )(S n - S 33 ) 
+B 1 (N 1 , N 2 )(S 13 + S 31 ) + B 4 (N U N 2 )S 22 , 

where B^Ni, N 2 ), i — 0, 1, 3, 4, are given by 



(111) 



B (N U N 2 ) = I( £l + £ 3 ) + E U(N t + |) - ^g^Ni) 



B 1 (N 1 ,N 2 ) 



i=i 



'' 1 +-t)9i(N 1 )g 2 {N 2 ), 



2VA1 1 A 2 - 



Bs(Ni, N 2 ) = \(e 1 - e 3 ) + § LZi - cu 2 Z 2 + 



B 4 (N U N 2 ) =e 2 +J2 u;^ + 

i=i L 1 



gl(Ni) gj(N 2 ) 



Ai 



A 2 



(112) 



where ft(^) = ft^) + 3, i = 1, 2. The Hamiltonian H" is a function of N x and 
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N 2 , and composed of two independent parts: H'{ = B (N 1 , N 2 )(S 11 + S 33 )+B 3 (N 1 , N 2 )(S U - 
S 33 ) + Bx(N 1: N 2 )(S 13 + S 31 ) is a two-level system containing a two-mode cavity; the other 
part is H'2 = B^Ni, N 2 )S 22 . The eigenkets of H" are directly written as 



\* + (ki,k 2 ))" 



^ COS f fc 2 ) ^ 



sin||fci, k 2 ) ^ 



^ - sin fc 2 ) ^ 



|*"(A;i,A; 2 ))" 



V 



o 



COS fl&i, fc 2 ) 



(113) 



/ 



where = 9(k 1 , k 2 ) = cos" 1 {^(fci, k 2 )/ [B^h, k 2 ) + flf^ k 2 )} 1/2 }. The eigenvalues cor- 
responding to these two kets are respectively 

1/2 



E ± (k 1 , k 2 ) = B (ki, k 2 ) ± Bf(h, k 2 ) + Bi(h, k 2 ) 



The third eigenket set of H" is determined by if 2 , 



(114) 



' ^ 



|* 3 (A;i,A; 2 ))" 



) 

the eigenvalues are simply 



V 



115) 



E 3 (h,k 2 ) = B 4 (h,k 2 ). 



(116) 



Using V\ we can turn these ket vectors to the original frame H' . \^ 3 (k±, k 2 ))" are invariant 
under V*. But ^(h, k 2 ))" turn out to be 



l^i,**))' 



^ cos + h, k 2 ) ^ 


y sin fc 2 + / 2 ) 



^ - sin + /i, fc 2 ) ^ 


y cos A; 2 + Z 2 ) 



(117) 



The results (114) and (117) indicate the existence of Rabi frequency 



v(fa, k 2 ) = Bl{k u k 2 ) + Bl{k u k 2 ) 



1/2 



(118) 



in the system H' . Therefore, some interesting effects appearing in the usual two-level system 
can also appear here. However, the above comment (ii) has indicated that these ket vectors 
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(117) are in fact the "dressed" ket vectors, induced by the unitary transformation Q as Eq. 
(98). For the viewpoint of experiment, it is meaningful to find these ket vectors in the "bare" 
form. Up to the first order approximation, we can use Eq. (105), and find that \^/ 3 (ki, k 2 ))" 
maintain the same form, but the vectors |^ r± (A;i, k 2 ))' turn out to be 

\ 



^ cos f |&4 + Zi, k 2 ) ^ 
Xi\h, k 2 ) 
sin k 2 + l 2 ) 



\y-(h,k 2 )) = 



^ - sin + h, k 2 ) 
X2\ki,k 2 ) 
cos f A; 2 + Z2) 



(119) 



where 



Xi 



pi 



2 ' 



(120) 



X2 = 



^V^sinf + ^V^cosf. 



The normalization factors are not included in the above ket vectors. Now we can conclude 
that only for the ket vectors as Eq. (119), the Rabi oscillation can appear in experiment. 

Recall that H' results from the first order approximation of H under the condition that 
(^aT^)' ( ^a^ ) ^ ^i' ^2- We now consider the following two special cases to reduce the 
above eigenvalues Eq. (114): 

First, for the case that is often used in the literature: 



Ai = A 2 = A, 

these E ± (ki, k 2 ) are reduced into 



E ± (k 1 ,k 2 )= 1 -(s 1 + e 3 )+j2 
1 i=i 



± 



9i(h) - al{k 2 ) 



2A 



The second case is contrary to Eq. (121) that 
we may accordingly reduce E ± (k 1 , k 2 ) into the following form 



E ± (k 1 ,k 2 ) = -(e 1 + e 3 ) + J2 

1 i=l L 



1 

±- 

2 



A A 1 9l{ki) , , 1 1^1(^1)^2(^2) 
A 2 - Ai + — + {— + — )- 



Ai 



•Ai A 2 / 2(A 2 -A 1 ) 



121) 



(122) 



(123) 



(124) 
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3.3.2. Exact solution 



We know that the approximate solution to H in the above subsection is based on the 
condition 

<«, « A„ A, (125) 

However, the quantities gi(Ni) = pi(Ni)[(Ni + Zj)!/7Vj!]5, 2 = 1, 2, increase with the photon 
numbers TVj for certain densities, which means that the condition (125) does not hold for 
the states of large photon numbers . In this subsection, we will solve the Hamiltonian (96) 
exactly. 

We employ again the right-unitary operator V in Eq. (106) as a transformation. It is 
easily checked that the Kernel of V given by Eq. (107) is still a set of the eigenkets of H. 
Since the Kernel is invariant under the unitary transformation Q [Eq. (98)], the eigenvalues 
of H corresponding to these eigenkets are still the same as those of H' [Eq. (100)]. Based 
on these results, we can use V as a unitary transformation to H to obtain its remaining 
subspace. Let 

H = VHV\ (126) 
a direct calculation gives 

H = H + Hi, (127) 
where 

3 2 

i7o = J>/3 + 5>i(Wi + lj/3), (128) 
i=i j=i 

and Hi is a traceless matrix, 
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/ 



h gi(Ni) o 

H l = g^N,) -(/i + / 2 ) 92{N 2 ) 
92(N 2 ) h 

where g\ (Ni) and g 2 (N 2 ) are given in Eq. (112), and 



V 



(129) 



/ 



h = |(2ei - £ 2 - e 3 + 2wiZi - cu 2 Z 2 ) = §(A 2 - 2Ax), 
h = |(2e 3 - £1 - £2 + 2w 2 / 2 - W1Z1) = |(Ai - 2A 2 ). 



(130) 



The Hamiltonian now becomes a function of photon numbers Aq and N 2 , where the creation 
and annihilation of the photon in the transitions between the atomic states are erased by the 
transformation V. Therefore, in the new frame, the transitions can only happen between the 
atomic states having the same photon number Aq (and N 2 ). For a ket vector with photon 
numbers k\ and k 2 



V 

fa 



(131) 



where fa = fa(k±, k 2 ), i = 1, 2, 3, i7 becomes the usual stationary three-level system, where 
the eigenvalues of H 1 are determined by the equation det \H 1 — A| =0, namely, 



A 3 + \p! + qi = 0, 
where 

Pl = - [^(fcx) + £ 2 (A; 2 ) + A 2 + /| + f\f 2 ] , 

. 5i = fi9l(h) + M 2 (A:i) + (/1 + / 2 )/i/ 2 . 
The solutions of this equation are 21 , 

Ai = Ai(fci,fc 2 )=C + + C_, 
1 A 2 , 3 = A 2 , 3 (Aq, k 2 ) = \{C + + C_) ± 4^(C + - C_), 



(132) 



(133) 



(134) 
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where C± = C±(k\, k 2 ) are given by 



2«i ± + i^i 



1/2 



(135) 



Combining the above \ with the term if in H, we obtain the exact eigenvalues of the 
three-level system 



£i(*4, fc 2 ) = E £i/3 + E ^-(fcj + Zj/3) + Ai, 
i=i j=i 

£ 2 ,3(£;i, fc 2 ) = E £i/3 + E Uj(kj + lj/3) + A 2 , 3 . 
i=i j=i 



(136) 



We assume that there is no degeneracy in the energy spectrum. Then, the eigenstates of 
H corresponding to the above eigenvalues are orthogonal mutually, which are given by 




1,2,3, 



(137) 



where ^ are the normalization factors: £j = 



1 + fe) 2 + te) : 



1/2 



. Using Theorem 



IV, we obtain the eigenkets in the original frame as: 



' '(A)l*i + Z i>**> N 



|*i(fci,A;2)> 



i = 1,2,3. 



(138) 



/ 



\ki, k 2 ) 
{(^)\ki,k 2 + l 2 ) 

Up to now, we have solved exactly the system of A-configuration atoms interacting with 
a two-mode cavity field by the RUT method. One may notice that the above procedure 
is almost the same as that in subsection 3.2, this means that to the RUT method, that 
the system contains one-mode cavity or two-mode cavity has no difference, even though the 
physical meanings of solutions are quite different. Using the expression of f\ and f 2 , one 
can expand the above solution according to the condition of large detunings, and compare 
the result with the approximate solution obtained in last subsection. Based on the above 
solution, one can further study the physical effects in this model. 
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In the literature, ones study the three-level system where the cavity fields change with 
time, i.e., the nonstationary case. We should point out that one can follow Theorem V to 
treat this nonstationary system, and the procedure is similar to that in above Subsection 
3.1. 

3.4. Application of RUT to the atoms with other configurations 

The above approach concentrates on the system where the atoms have A-configuration en- 
ergy levels only. Recently, ^-configuration atoms were found to exhibit some new interesting 
effects 24 . In principle, under the condition that the detunings are very large, a system of 
^-configuration atoms interacting with cavity fields can be treated by the same way as the 
above approximate treatment to A-configuration atoms. Unfortunately, for =-configuration 
atomic system, this treatment becomes invalid. 

However, as we have pointed out that the detunings are not always very large for various 
atomic systems. Moreover, gi(Ni) in Eq. (112) always increase with the photon number N iy 
but the detunings are invariant with JVj. Therefore, the perturbation treatment introduced 
above is valid only to the situation of large detunings and low photon number state. One may 
notice that the exact treatment of three-level atomic system by RUT method is evidently 
simpler than the perturbation treatment, and suitable to various situations. In this subsec- 
tion, we briefly show how to apply RUT method to the systems of V and =-configuration 
atoms. 

We first look at a system of U-configuration atoms interacting with one-mode cavity field 
shown in Fig. 2.b. The generalized Hamiltonian is 

3 

H = J2 £ i S a + + [pi(N)a h S 12 + a^p 2 (N)S 2 3 + /i.e.], (139) 
where we assume l\ < l 2 . To solve this Hamiltonian, we here construct an operator U 
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u 



/ F h _ h 0^ 
F h 
1 



(140) 



where F^-^ and Fi 2 are two phaser operators. U belongs to RUT, and the Kernel of U is 



K = i 



l#(*i)> 





v I 



, l*S(* 2 )> 



I o ^ 
V / 



, ki<l 2 — h, k 2 < l 2 



(141) 



We can take |\&) = XilV'i(^i)) + X2\^ 2 {k 2 )) to check whether |\&) is the eigenket of H. With- 
out difficulty, we find the following sets of vectors 





cos \ki) 




— sin \ki) 




( \ 


liM*0> = 


sin%|A;i + Zi) 




cos%|A;i + Zi) 


, \Uk2)) = 






V ° J 









{ J 



, (142) 



are the eigenkets of H, where k± < l 2 — li, k 2 < Z 1; 9 kl = tan -1 [2^ (ki)/(ei — e 2 — uli)], 
and g^iki) = pi(k 1 ) [(ki + li)\/k^] 1 ^ 2 . The eigenvalues corresponding to above eigenkets are 
obtained to be 



E^h) = \{e l + e 2 + uh) + uh + ^\(e l -e 2 -ujl l ) 2 + g h (k l )\ 
E^h) = i(ei + e 2 + wZi) + wfci - ^(ei-ea-^iP + ^Cfci) 2 , 
£3(^2) = e 2 + uk 2 . 



(143) 



On the other hand, one can check that each element in the Kernel can be written as a linear 
combination of eigenkets in Eq. (141). These results indicate that U can be applied to H 
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to obtain the its remaining subspace. A simple calculation gives 



H' = UHU ] = 



( e 1 + cu(N + l 2 -l 1 ), g h (N + l 2 -h), 



g h (N + l 2 -h), e 2 + u(N + l 2 ), g h (N) 



V 



0, 



g l2 (N), e 3 + ujN) 



(144) 



The Hamiltonian now becomes a matrix function of photon number N. One can follow the 
procedure in the above subsections to obtain the eigenvalues and eigenkets of H'. We omit 
these here. 

We now look at a system of three-level atoms with =-configuration energy level, which 
interact with a one-mode cavity as shown in Fig. 2.c. The Hamiltonian is 

3 



H = Y,£iSu + uja^a + [ Pl {N)a h S l2 + p 2 {N)a h S 23 + h.c], 



(145) 



i=l 



To solve the Hamiltonian, we introduce an operator matrix V as 



^10 N 



V 



F h 
v F h+ i 2 

V evidently belongs to RUT. The Kernel of V is 



(146) 





( ] 




( \ 


> 


< l$(*l)> = 


l*l> 


, l«fe)) = 





, ki < h, k 2 < h + l 2 > 




V 1 




V ^ ) 





Using the same method as in the above case of U-configuration, we obtain that this Kernel 
is isomorphic with to subset of the eigenkets of H, where the eigenkets are 
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liM*i)> 



cos -2~\ki) 
sin + l 2 ) t 



, |^ 2 (fci)> 



V 



-sin-^lfci) 
cos%|A:i + / 2 ) 



, |^i(fc 2 )> 



/ 



' ^ 







, (148) 



where k\ < h, k 2 < l 2 , and 9 kl = tan 1 \2gi 2 (ki) / (e 2 — e 3 — ool 2 )}. The eigenvalues corre- 
sponding to above eigenkets are 



E i(ki) = i (£2 + ^3 + wh) + uki + ^\{e 2 -ez-uol 2 y + g l2 {k 1 )\ 



E^h) = |(e 2 + e 3 + + uh - ^Jl(e 2 - e 3 - ul 2 ) 2 + g^h) 2 , 
El(k 2 ) =e 3 + uk 2 . 



(149) 



One can prove that an arbitrary element in the Kernel K can be expressed by the above 
\ijji(k)), \ip 2 (k)) and \ifj 3 (k)} in linear form. This means that K is isomorphic with one set 
of the eigenkets of H . Therefore, the operator U satisfies Theorem IV, and can be applied 
to the remaining eigenket set of H. We obtain 



H' = VHV ] 



\ 



e.+uN, g h (N), 



g h (N), e 2 + u(N + h), g h (N + h) 



\ 0, g l2 (N + h), e 3 + u;(N + h + l 2 ) J 



(150) 



H' is a matrix, whose matrix elements are functions of photon number N. H' can be solved 
by the regular way. We omit these here. 

The above approach shows that all the A, V and =-configuration atoms interacting with 
quantized cavity fields can be treated unitedly by RUT method, where the right-unitary 
operators are simply diagonal matrices constructed by phaser operators. 

In the last part of this section, we would like to show how to apply the RUT method 
to the atom-radiation interaction system where the atoms have higher (> 3) energy level. 
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For an example, we here choose a simple case: four-level atoms with configuration shown in 
Fig. 3. Under the rotating- wave approximation, the generalized Hamiltonian of the atoms 



interacting with one-mode cavity field is 

4 

H = + ujata + [pia h S 12 + p 2 a h S 23 + p 3 a h S u + h.c], 



(151) 



i=l 



for simplicity, we assume pi i — 1, 2, 3 be constants. To solve this Hamiltonian, we construct 
a right unitary operator as 
^10 
F h 
F h+l2 
^ F h+h+h J 

The Kernel of Va is 



V 4 



(152) 



K = I 





( \ 






< l^?(*l)> = 













v , 



' o N 



V / 



l^3°(*3)> 



' ^ 



V ^ ) 



> . 



(153) 



where k\ < h, k 2 < h + l 2 , k 3 < l\ + l 2 + l 3 . Within these ket vectors, the four-level 
system is reduced into a three-level case, one can further check out that these ket vectors 
are isomorphic with a subset of the eigenkets of H. Which means that we can use Va as a 
unitary operator to the other subspace of H. Let H' = V4HV4, then 



H' 



e 1 + uN, 9l (N), 0, 
9l (N), £2 + u(N + h), g 2 (N) 

0, g 2 {N), e 3 + u(N + h + l 2 ) 

0, 0, £3 (AO, 
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^ 



S3 (AO 

Ea + uj{N + h + l 2 + / 3 ) J 



(154) 



where 9l (N) = Pl [(N + h)\ / N\]$ , g 2 (N) = p 2 [(N + h + l 2 )\ / \N + h)^ , g 3 (N) = p 2 [(N + 
h + h + ls)\/(N + li + hy]?- For the ket with a fix photon number N, H' is simply a 4 x 4 
constant matrix. Thus, its eigenkets and eigenvalues are easily obtained. 

One can further follow this method to solve other configurations of the four-level systems. 
These are omitted here. 



IV. Conclusion 

In conclusion, we have developed the right-unitary transformation (RUT) theory, and ini- 
tially discussed its applications in physics. The first part of this paper attributes to the 
theory. We found that the state space of any operator U e RUT (precisely W) is composed 
of two independent part: { |^°) } and { }, where { |\1/ ) } is called the Kernel of U, 
which satisfies U{ \^°) } = 0. On the other hand, in the subspace { l^ 1 ) }, U acts as a 
unitary operator. The properties of RUT such as semigroup, Kernel, etc., were discussed 
through several theorems. Based on these properties, we concluded a general way on how 
to apply the RUT to a physical system. For a physical quantity such as the Hamiltonian H, 
suppose its eigenstates as S = { , % — 1, 2, • • • oo }, then 

oo 

H = (155) 
i=i 

Ei are eigenvalues. If a subset of S: Si — { |<E»j) , % — 1, 2, • • • s } is checked to be isomorphic 
with the Kernel of U, then the supplement set of Si is evidently isomorphic with the unitary 
subspace of U . Therefore, 

oo 

H' = UHtf = ]T E^i&l (156) 

i=s+l 

where |<3>-) = U\&i). The new state space S' = { |$-) , % — s + 1, • • • oo } is proved to be 
complete. Equation (156) shows that the new frame H' has the same spectrum (precisely a 
subset of the spectrum) as H . 
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Based on the above results, in the second part of this paper we used the RUT method to 
deal with the systems of many-level atoms interacting with the quantized radiation fields, 
where the RUT are the matrices constructed by the well-known phase operators. We have 
studied two typical systems: One is about the Jaynes-Cummings models, which were found 
to exhibit some interesting effects, and have obtained much study in recent years. We solved 
a nonstationary generalized JC model, and found that atomic inversion of the system can 
be controlled through the external source. Another system carefully studied is about the 
interaction of the three-level atoms with one or two-mode cavity field. This system has 
been widely applied in various contexts of quantum optics such as lasing without inversion, 
electromagnetically induced transparency, etc. This paper provides a unified method for 
these topics. 

We would like to point out that the RUT method can be applied to some simplified 
quantum electrodynamic (QED) system, such as the photon-electron, and phonon-electron 
interactions. It will be discussed in the forthcoming presentation. 
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Appendix 



For the initial state Eq. (63), 



z 



k 




(1) 
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Then, 



From Eq. (51), we have 

uj+ + ujk + f - 9l (k) 



9i(k) uj^+ujk + f : 

and 

out + uok + 



2 _ ^fc I u/.v | 2 
LU^ +LUk + ^f 



Using these results, we rearrange af as 

a k = [ , ,- ; ,+ /(*)> «*=- ,.- - 

with these expressions, we directly obtain Eqs. (64) and (65) . 
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FIGURES 

Fig. L (a) shows that in the frame H, the transition from the atomic state 2 to state 1 is achieved 
by creating photons, (b) shows that in the transformed frame H' , there is no photon created or 
annihilated in the transition from 2 to 1. 

Fig. 2. (a) represents the transitions of the atoms with A configuration energy levels, (b) indicates 
the V configuration, and (c) the = configuration. 

Fig. 3. Four-level atoms with ladder-configuration energies. 
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